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ABSTRACT 

o 

O ■ There is now good evidence linking short-hard GRBs with both elliptical and 

spiral galaxies at relatively low redshifts, redshift of about 0.2. This contrasts 
with the average redshift of about 2 of long-duration events, which also occur 
only in star-forming galaxies. The diversity of hosts is reminiscent of type la 
supernovae, which are widely (but not universally) believed to originate from 
the coalescence of white dwarfs. By analogy, it has been postulated that short- 
hard bursts originate from neutron star mergers. Mergers, as well as stellar 
core-collapse events (type II SNe and long-duration GRBs) are accompanied by 
long-lived sub-relativistic components powered by radioactive decay of unsta- 
ble elements produced in the explosion. It is therefore interesting to explore 
whether short duration events also have ejecta powered by radioactivity (i.e. 
that are supernova-like). Observations already inform us that any supernova 
like component in the first few well studied short hard bursts must be fainter 
than those typical of type la or core-collapse supernovae. Rather than refer to 
weaker supernova-like component as "mini-super nova", an etymologically in- 
defensible term, I use the term macronova. I investigate the observability of 
j> ! macronovae powered by neutron decay and by radioactive Nickel. Separately, I 



note that a macronova will reprocess energetic emission arising from a long lived 
central source. I find that surprisingly interesting limits on the basic parameters 
of macronovae can be obtained provided observations are obtained with current 
10-m class telescopes over a range of one hour to one day following the burst. 
Detection or even strong constraints provided by macronova observations can be 
expected to yield the most direct clue (barring data from futuristic gravitational 
wave interferemeters) to the origin of short hard bursts. It is this prize that, over 
the course of coming few years, will undoubtedly motivate astronomers to mount 
ambitious macronova campaigns. 

Subject headings: gamma rays:bursts — supernovae: general — stars:novae 
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1. Background 

There is currently a great deal of excitement and sense of great progress in our under- 
standing of the mysterious class of short duration, hard spectrum gamma-ray bursts. GRB 
050509b was the first short hard burst to be localized to moderate precision (< 10"). This 
was made possible thanks to the prompt slewing of Swift and subsequent discovery of a 
rapidly fading X-ray source (Gehrels et al. 2005). The X-ray afterglow was noted to coincide 
with a nearby cluster of galaxies (Gehrels et al. 2005) and lay in the outskirts of a bright 
elliptical galaxy at z = 0.22 (Bloom et al. 2005). Fortunately the next two short hard bursts 
were localized to arcsecond accuracy thanks to the discovery of their X-ray, optical or radio 
afterglows and, as a result, greatly advanced the field of short hard bursts. GRB 050709 was 
detected by the High Energy Transient Explorer (Villasenor et al. 2005). The optical (Hjorth 
et al. 2005) and X-ray (Fox et al. 2005) afterglow places this burst on a star-forming galaxy 
at z = 0.16 (Fox et al. 2005). GRB 050724 was detected by Swift (Covino et al. 2005). The 
radio and optical afterglow localization places this event on an elliptical galaxy (Berger et al. 
2005) at z = 0.26 (Prochaska et al. 2005). 

Collectively these three short bursts indicate both elliptical and spiral host galaxies. 
These associations lend credence to other less well localized associations e.g. Nakar et al. 
(2005). It is however a recent development that is of potentially greater interest to this 
paper, namely the claim by Tanvir et al. (2005) of a population of short bursts that are even 
closer. 

It is intriguing that the hosts of short hard bursts are as diverse as those of la supernovae. 
A model (popular in the past but not ruled out by observations) that is invoked for la 
explosions involves double degenerate binary systems which at some point coalesce and then 
explode. An early and enduring model for short hard bursts is the coalescence of two neutron 
stars (or a neutron star and a black hole). The coalescence is expected to produce a burst 
of neutrinos, of gamma-ray and also eject neutron- rich ejecta (Eichler et al. 1989). 

Several years ago Li & Paczyhski (1998) speculated the ejecta would result in a supernova- 
like explosion i.e. a sub-relativistic explosion with radioactivity. The purpose of this paper 
is to revisit this topic now that the rough distance scale has been determined. Furthermore, 
the possibility that some events may be even very close makes it doubly attractive to develop 
quantitative models of associated supernova-like explosions. 

Coalescence models are not silent on the issue of associated non-relativistic outflows (e.g. 
Janka & Ruffert 2002; Rosswog 2005). Indeed, there appears to be a multiplicity of reasons 
for such outflows: tidal tails (which appear inevitable in all numerical models investigated to 
date), a wind driven by neutrino emission of a central massive neutron star and explosion of 
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a striped neutron star. To this list I add another possible mechanism: ejection of the outer 
regions of the accretion disk owing to conservation of angular momentum. 

The composition of the ejecta is a matter of fundamental theoretical interest but also of 
great observational consequence. For the explosion to be observationally detected the ejecta 
must contain a long-lived source of power. Otherwise the ejecta cools rapidly and is not 
detectable with the sensitivity of current telescopes. Indeed, the same argument applies to 
ordinary SN. Unfortunately, theoretical models of coalescence offer no clear guidance on the 
composition of the ejecta although the current prejudice is in favor of neutron rich ejecta 
(see Freiburghaus et al. 1999). 

The paper addresses two potentially novel ideas: an explosion in which the decay of 
free neutrons provides the source of a long lived source of power and the reprocessing of 
the luminosity of a long lived central source into longer wavelength emission by the slow 
moving ejecta. A less speculative aspect of the paper is that I take a critical look about 
photon production and photon-matter equilibrium. These two issues are not important for 
SN models (and hence have not been addressed in the SN literature) but could be of some 
importance for lower luminosity and smaller ejecta mass explosions. Along the same vein, 
I have also investigated the transparency of the ejecta to 7-rays - an issue which is critical 
given the expected low mass and high speed of the ejecta (again relative to SN). 

The two ideas discussed above (a neutron powered MN and long lived central source) 
are clearly speculative. However, the model presented here include the essential physics of 
such explosions but are adequate to explore the feasibility of detecting associated explosions 
over a wide range of conditions. 

Now a word on terminology. It is clear from the observations of the three bursts and their 
afterglows that any accompanying sub-relativistic explosion laced with or without radioactive 
isotopes is considerably dimmer than typical supernovae (la or otherwise; Bloom et al. 
2005; Hjorth et al. 2005; Fox et al. 2005; Berger et al. 2005) The word "mini supernova" 
may naturally come to one's mind as an apt description of such low luminosity explosions. 
However, the juxtaposition of "mini" and "super" is not etymologically defensible, and will 
only burden our field with more puzzling jargon. After seeking alternative names I settled 
on the word macronova (MN) 2 - an explosion with energies between those of a nova and a 
supernova and observationally distinguished by being brighter than a typical nova (My ~ 
—8 mag) but fainter than a typical supernova (My ~ —19 mag). 



2 This word was suggested by P. A. Price. 
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2. The Physical Model 

All short hard burst models must be able to account for the burst of gamma-ray emission. 
This requires ultra-relativistic ejecta (Goodman 1986; Paczynski 1986). Here, we are focussed 
entirely on possible ejecta but at sub-relativistic velocities. The fundamental parameters of 
any associated macronova is the initial internal energy (heat), E , the mass (M cj ) and the 
composition of the sub-relativistic ejecta. Given that the progenitors of short hard bursts are 
expected to be compact the precise value of the initial radius, Ro, should not matter to the 
relatively late time (tens of minutes to days) epochs of interest to this paper. Accordingly, 
rather arbitrarily, Rq has been set to 10 7 cm. There is little guidance on E and M e j but 
M ej ~ 1O" 4 M to 1(T 2 M (and perhaps even 0.1 M & ) have been indicated by numerical 
studies (Janka & Ruffert 2002; Rosswog 2005). Based on analogy with long duration GRBs, 
a reasonable value for E is the isotropic 7-ray energy release of the burst itself. We set the 
fiducial value for E to be 10 49 erg. 

It appears to me that there are three interesting choices for the composition of the ejecta: 
a neutron rich ejecta in which the elements decay rapidly (seconds; Freiburghaus et al. 1999), 
an ejecta dominated by free neutrons and an ejecta dominated by 56 Ni. Isotopes which decay 
too rapidly (e.g. neutron rich ejecta) or those which decay on timescales much longer than 
a few days will not significantly make increase the brightness of the resulting macronova. 
A neutron- and an 56 Ni-MN are interesting in that if such events exist then they are well 
suited to the timescales that are within reach of current observations, 25th magnitude 3 on 
timescales of hours to days. 

The explosion energy, E , is composed of heat (internal energy) and kinetic energy of the 
ejecta. The heat further drives the expansion and so over time the internal energy is converted 
to additional kinetic energy. We will make the following, admittedly artificial, assumption: 
the initial heat is much larger than the initial kinetic energy. The final kinetic energy is 
thus E = l/2M e jVg where v s is the final velocity of ejecta. Consistent with the simplicity 
of the model we will assume that the expanding ejecta is homogeneous. The treatment in 
this paper is entirely Newtonian (unless stated otherwise) but it is more convenient to use 
f3 s = v s /c rather than v s . 

At any given instant, the total internal energy of the expanding ejecta, E, is composed 



3 The following note may be of help to more theoretically oriented readers. The model light curves 
presented here are for the Johnson / band (mean wavelength of 0.8 /i) / = corresponds to 2550 Jy. The AB 
system, an alternative system, is quite popular. This system is defined by m(AB) = —2.5 log 10(/„) — 48.6 
where /„ is the flux in the usual CGS units, erg cm~ 2 s _1 . This corresponds to a zero point of 3630 Jy at 
all frequencies. Thus 25th magnitude is a few tenths of /xJy. 
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of a thermal term, E t ^ arising from the random motion of the electrons (density, n e ) and 
ions (density, rii) and the energy in photons, E ph : 

E/V = -rn(Z + l)kT + aT\ (1) 

Here, V = 4nR 3 /3, N t = M e j/(iro 5 ), = Ni/V and n e = Zrii. For future reference, the 
total number of particles is N — N^Z + 1). Implicit in Equation 1 is the assumption that 
the electron, ion and photon temperatures are the same. This issue is considered in some 
detail in §4. 

The store of heat has gains and losses described by 

E = e(t)M ej - L{t) - 4nR(t) 2 Pv(t) (2) 

where L(t) is the luminosity radiated at the surface and e(t) is heating rate per gram from 
any source of energy (e.g. radioactivity or a long-lived central source). P is the total pressure 
and is given by the sum of gas and photon pressure: 

P = m(Z + l)kT + aT 4 /3. (3) 



As explained earlier, the ejecta gain speed rapidly from expansion (the 4nR 2 Pv work 
term). Thus, following the initial acceleration phase, the radius can be expected to increase 
linearly with time: 

R(t) =R + fact. (4) 

With this (reasonable) assumption of coasting we avoid solving the momentum equation and 
thus set v = v s in Equation 2. 

Next, we resort to the so-called "diffusion" approximation (see Arnett 1996; Padman- 
abhan 2000, volume II, §4.8), 

L = E ph /t d , (5) 

where 

t d = BkM c] /cR (6) 

is the timescale for a typical photon to diffuse from the center to the surface. The pre- 
factor B in Equation 6 depends on the geometry and, following Padmanabhan {ibid), we set 
B = 0.07. k is the mass opacity. 

The composition of the ejecta determines the heating function, e(t), and the emission 
spectrum. For neutrons, the spectrum is entirely due to hydrogen. Initially the photosphere 
is simply electrons and protons and the main source of opacity is due to Thompson scattering 
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by the electrons. The mass opacity is k — a T /m H = 0.4cm 2 g _1 ; here m# is the mass of a 
hydrogen atom and Ot is the Thompson cross section. Recombination of protons and elec- 
trons begins in earnest when the temperature reaches 20,000 K and is complete by 5,000 K. 
With recombination, the opacity from electrons disappears. Based on models of hydrogen 
rich SNe we assume that the critical temperature, the temperature at which electron opacity 
disappears, is Trc = 10 4 K. 

In contrast, for Nickel (or any other Z ^> 1 element) ejecta, the spectrum will be 
dominated by strong metal lines (like la supernovae) with strong absorption blueward of 
4000 A. Next, Z/A ~ 0.5 and k = 0.2 cm 2 g" 1 . In this case, based on models for la SNe 
spectra (Arnett 1996), I assume that the photosphere is entirely dominated by electrons for 
T > T RC = 5 x 10 3 K. 



The Thompson optical depth is given by 

ZMci/mH Ra T = 20(MJ10- 3 M Q )(Z/A)R U "; (7) 



es A 47ri? 3 /3 

here, we use the notation Q x as short hand notation for a physical quantity normalized to 
10 s in CGS units. Thus the ejecta, for our fiducial value of M e j, remain optically thick 
until the size reaches about 10 14 cm. However, as noted above, electron scattering ceases for 
T < Trc following which the use of Equation 7 is erroneous. The reader is reminded of this 
limitation in later sections where the model light curves are presented. 

With the great simplification made possible by Equations 4 and 5, the RHS of Equation 2 
is now a function of E and t. Thus we have an ordinary differential equation in E. The 
integration of Equation 2 is considerably simplified (and speeded up) by casting P in terms of 
E (this avoids solving a quartic function for T at every integrator step). From Equation 3 we 
find that photon pressure triumphs over gas pressure when the product of energy and radius 
ER> x = 5 x 10 55 ((Z + l)/A) 4 / 3 (M cj /10- 2 M ) 4 / 3 ergcm. The formula P = (1 + x/(x + 
ER))E/3V allows for a smooth transition from the photon rich to photon poor regime. 

Applying the MATLAB ordinary differential equation solver, odel5s, to Equation 2 I 
obtained Bona logarithmic grid of time, starting from 10 ms to 10 7 s. With the run of E 
(and R) determined, I solved for T by providing the minimum of the photon temperature, 
(E/aV) 1 / 4 and the gas temperature, 2E/(3iV/c B ) as the initial guess value for the routine 
f zero as applied to Equation 1. With T and R in hand E pil is easily calculated and thence 
L. 

The effective temperature of the surface emission was computed using the Stefan- 
Boltzman formula, L = ATTR 2 aT e s A . The spectrum of the emitted radiation was assumed 
to be a black body with T eff . Again this is a simplification in that Comptonization in the 
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photosphere has been ignored. 



3. Pure Explosion 



In this section we consider a pure explosion i.e. no subsequent heating, e(t) = 0. If 
photon pressure dominates then P = 1/3(E/V) and an analytical formula for L(i) can be 
obtained (Padmanabhan, op cit; Arnett 1996): 



here, th = Ro/v s is the initial hydrodynamic scale, td(0) = B(kM^/cRq) is the initial diffu- 
sion timescale and L = i? /^(0). However, for the range of physical parameters discussed 
in this paper, gas pressure could dominate over photon pressure. Bearing this in mind, 
I integrated Equation 2 but with P = 2/3(E/V) and found an equally simple analytical 
solution: 



The analytical formula allow us to get some insight into the overall behavior of the 
luminosity evolution. First, we note that th = Ro/v s = OS/S' 1 ms is much smaller than 
td(0) = 6.2 x 10 3 (M C j/10~ 3 M Q ) yr. The internal energy, E, decreases on the initial hydrody- 
namical scale and immediately justifies our coasting approximation. Next, the duration of 
the signal is given by the geometric mean of td(0) and th and is oc (M c Jv s ) l l 2 but indepen- 
dent of Rq. The duration is not all that short, ~ 0.3(M cj /l(r 3 M & ) l / 2 ((3 S /Q.l)- 1/2 d. Third, 
the peak emission, Eo/td(0) = (3 2 c 3 Rq/(2Bk), is independent of the mass of the ejecta but 
directly proportional to Rq and the square of the final coasting speed, v 2 . 

Unfortunately, as can be seen from Figure 1 (bottom panel), model light curves 4 for a 
macronova located at z — 0.2 is beyond current capabilities in any band even in the best 
case (high (5 S ). This pessimistic conclusion is a direct result of small Rq, small M c - S and great 
distance (~ 1 Gpc) for short hard bursts. 



4 Thc model light curves presented here, unless stated otherwise, are for a luminosity distance of 0.97 
Gpc which corresponds to z = 0.2 according to the currently popular cosmology (flat Universe, Hubble's 
constant of 71 kms -1 Mpc -1 ). The two wavebands discussed here are the Optical (corresponding to restframc 
wavelength of 8140 A/(l + z)) and the UV band (corresponding to rcstframe wavelength of 1800 A/(l + z)). 
These two bands are chosen to represent ground-based I-band (a fairly popular band amongst observers 
given its lower immunity to lunar phase) and one of the UV bands on the Swift UV-Optical telescope. The 
time axis in all figures has not been stretched by 1 + z. 
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Fig. 1. — (Top) The luminosity, L(t), for an explosion with no heating source. The model pa- 
rameters (M e j and (3 S ) are indicated on the title line. The thick line is obtained from numerical 
integration (Equation 5) whereas the dotted line is the (photon dominated) analytical formula 
(Equation 8). The numerical curve tracks the analytical curve (apart from a scaling of 0.6); the 
two disagree as the MN evolves (see text). (Bottom) The optical and UV broad-band fluxes (vf v ) 
expected for a macronova located at a redshift of z = 0.2. The dotted line is the bolometric flux, 
L(t)/(4vrL> 2 ) where D = 0.97 Gpc is the distance. 



The situation is worse for lower shock speeds. A lower shock speed means a lower pho- 
ton temperature and thus lower photon pressure. The luminosity decreases by an additional 
factor oc (t/t^)' 1 (Equation 9). Next, the internal energy is shared equitably between pho- 
tons and particles and a lower photon density means that a larger fraction of the internal 
energy is taken up the particles. Indeed, one can see a hint of the latter process in Figure 1 
where the numerical curve (proceeding from start of the burst to later times) is increasingly 
smaller than the analytical curve. This is entirely a result of equipartition of E between 
photons and particles and this equipartition is not accounted for by in deriving Equation 8. 



4. Heating by Neutron Decay 

On a timescale (half lifetime) of 10.4 minutes, a free neutron undergoes beta decay yield- 
ing a proton, a mildly relativistic electron with mean energy, 0.3 Me V) and an antineutrino. 
Thus the heating rate is entirely due to the newly minted electron (see Appendix) and is 



e n (t) = 3.2 x lO^ergg^s -1 . 



(10) 
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Even though there are two unconstrained model physical parameters, M ej and (3 S , a 
constraint may be reasonably imposed by the prejudice that the macronova energy is, at 
best, comparable to the gamma-ray energy release which we take to be 10 49 erg (Fox et al. 
2005). Within this overall constraint I consider two extreme cases for the neutron ejecta 
model: a high velocity case, f3 s = 0.3 and a low(est) velocity case, f3 s = 0.05. The mass of 
the ejecta was chosen so that E ~ 10 49 erg. The range in (3 is nicely bracketed by these two 
cases. The escape velocity of a neutron star is 0.3c. The energy released per neutron decay 
of 0.3 MeV is equivalent to (3 S = 0.025. Clearly, final ejecta speeds below this value seem 
implausible (even such a low value is implausible especially when considering that the ejecta 
must first escape from the deep clutches of a compact object). By coincidence, as discussed 
below, these two cases nicely correspond to a photon dominated and gas pressure dominated 
case, respectively. 



Neutron Decay: Mej= 3.0e-04 0=0.300 



Neutron Decay: Mej=1.0e-02 
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Fig. 2. — (Top) Internal energy (E) with (solid line) and without (dashed line) neutron decay 
heating. (Bottom) The interior temperature obtained by solving Equation 1 (solid line) given E 
and radius, R. The dash-dotted line is the temperature one would obtain if all the internal energy 
was in particles whereas the dotted line is that obtained if the photons dominated the internal 
energy. Each vertical pair of panels refers to a choice of model parameters: coasting speeds of 
P s = 0.3 (left) and (3 S = 0.05 (right). The explosion energy, E = l/2v 2 s M c] ~ 2 x 10 49 erg in both 
cases. "RC" marks the epoch at which the surface temperature falls below T c = 10 4 K and r = 10 
marks the epoch at which electron scattering optical depth is 10. "Ex" marks the epoch at which 
all the initial photons are radiated away. 



The decay of neutrons extends the hot phase of the fireball expansion (Figure 2). This is 
nicely demonstrated in Figure 2 (right set of panels) where we see that the decay of neutrons 
reverses the decrease in the internal energy. Indeed for the f3 s = 0.05 case neutron heating 
results in the pressure switching from being dominated by gas to a photon dominated ejecta. 
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It is this gradual heating that makes a neutron MN potentially detectable. 

I now carry out a number of self-consistency checks. To start with, implicit in Equa- 
tions 1 and 3 is the assumption that the ions, electrons and photons have the same temper- 
ature. From the Appendix we see that the electron-electron and electron-ion timescales are 
very short and we can safely assume that the electrons and ions are in thermal equilibrium 
with respect to each other. 

Next, I consider the time for the electrons to equilibrate with the photons. There are two 
parts to this issue. First, the slowing down of the ejected beta particle (electron). Second, 
the equilibration of thermal electrons with photons. The energetic electron can be slowed 
down by interaction with thermal electrons and also by interaction with photons. In the 
Appendix we show that the slowing down and thermalization timescales have essentially the 
same value: 



Thus when the interior temperature falls below (say) 2.7 x 10 4 K the photon-electron equili- 
bration time becomes significant, about 10 4 s. (However, by this time, most of the neutrons 
would have decayed and all that realistically matters is the photon-matter thermalization 
timescale). 

An entirely independent concern is the availability of photons for radiation at the surface. 
At the start of the explosion, by assumption, all the internal energy, E , is in photons and thus 
the initial number of photons is iV" p h(0) = E /(2.7kBTo) where E = VqoTq, Vq = 4tt/3RI 
is the initial volume and a is the Stefan-Boltzmann radiation constant (see Appendix). 
Photons are continually lost at the surface. Electron scattering does not change the number 
of photons. Thus an important self-consistency check is whether the number of radiated 
photons is consistent with the number of initial photons. However, as can be seen from 
Figure 2 the number of radiated photons exceeds the number of initial photons in less than 
ten minutes after the explosion (the epoch is marked "Ex"). Relative to the observations 
with large ground based telescopes this is a short timescale and so it is imperative that we 
consider the generation of new photons. 



For a hot plasma with A — Z — 1, the free-free process is the dominant source of 
photon creation. The free-free process is flat, f u oc exp(—hv/k B T) whereas the blackbody 
spectrum exhibits the well known Planck spectrum. Thus the free-free process will first 
start to populate the low-frequency end of the spectrum. Once the energy density, at a given 




(11) 



4.1. Free-free emission as a source of new photons 
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frequency, reaches the energy density of the black body spectrum at the same frequency, 
then free-free absorption will suppress further production of free-free photons. 



Mej= 3.0e-04 fS=0.300 



Mej=1.0e-02 |5=0.O50 




Fig. 3. — Photon diffusion (t^) timescale and photon-electron equilibration (i e , 7 ) timescales, rela- 
tive to the expansion time (i), as a function of t. The epoch at which these timescales match the 
expansion time are marked by an open square or open circle. 

Provided the free-free optical depth, tr(v) = Ra^iv) (Appendix) is ^ 1 for u ^ kT/h a 
simple estimate for the timescale over which the free-free process can build a photon density 
equal to that of the black body radiation field (of the same temperature) is given by 



ts 



aT 4 



e ff (n e ,nj,T) 



(12) 



here, eff(n e ,nj,T) is the frequency integrated free-free volume emissivity (see Appendix). 
However, electron scattering increases the effective optical depth of the free-free process 
(Rybicki & Lightman 1979, p. 33) 



t*{v) = \/ t s {v){t s {v) + r cs ) . 



(13) 



takes into account that the relevant optical depth for any emission process is the distance 
between the creation of a photon and the absorption of the photon. Electron scattering 
increases the probability of a newly minted photon to stay within the interior and thereby 
increases the free-free optical depth. 

Even after free-free process stops producing photons, the photons in the interior are 
radiated on the photon diffusion timescale (Equation 6). These two processes, the increase 
in the effective optical depth and the photon diffusion timescale, prolong the timescale of 
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Frgg-free timescale & Optical Depth Mej= 3.0e-04 beta=0.300 Free-free timescale & Optical Depth Mej= 1.0e-02 beta=0.050 
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Fig. 4. — Run of the free-free timescale, tg (top panel) and the effective free-free optical depth 
(r*) evaluated at normalized frequency x = hv/ksT = 2.7 (bottom panel). The timescales are 
normalized by time past the explosion. The free-free timescale is only meaningful when the effective 
optical depth is above unity (we have chosen the criterion) and the range of epochs for 

which this is the case is marked by a solid line. The epoch at which free-free emission is unable 
to keep the interior stocked with the blackbody photon density is marked by an open square 
(top panel). Right. The curves are for (3 S = 0.05 and M e j = 1O _2 M . The macronova never 
gets optically thin and hence the absence of squares. Left. The curves are for (3 S = 0.3 and 
M e j = 3 x 10~ 4 Mq. Seemingly td rises at late times but this artifact is of little consequence since 
there is no emission beyond 10 4 s. 

the macronova signal, making the macronova signal detectable by the largest ground-based 
optical telescopes (which unlike robotic telescopes respond on timescales of hours or worse). 

In Figures 3 and 4, I present various timescales (td, t-y, e , ts) and the frequency at which 
the effective free-free optical depth (Equation 13) is unity. For f3 s = 0.05, tg is always smaller 
than t and this means that the free-free process keeps the interior well stocked with photons. 
The remaining timescales are longer, approaching a day. 

For p s = 0.3 we find that r* as evaluated at x — hu/ksT = 2.7 falls below 2 at about 
2,000 s. Thus, for epochs smaller than 2,000 s we can legitimately evaluate tg. We find that 
ts exceeds t at about 1,000s (marked by "ff"). We consider this epoch to mark the epoch 
at which free-free emission "freezes" out in the sense that there is no significant production 
of free-free emission beyond this epoch. However, the photons in the interior leak out on 
the diffusion timescale which becomes comparable to the expansion time at t ~ 10 4 s. In 
summary, there is no shortage of photons for surface radiation for ejecta velocity as high as 
P s = 0.3 but only for epochs earlier than 10 4 s. 
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Fig. 5. — Optical and UV light curves for a macronova located at a redshift, z = 0.2. The physical 
parameters, M e j and (3 are displayed in the Figure. The symbols at the bottom of the figure are 
as follows: "Ex" (exhaustion of initial photons), "ff" (epoch beyond which free- free emission can 
no longer keep the photon energy density at the black body energy density), "7" (epoch at which 
electrons and photons decouple), "r es " (the electron scattering optical depth is 10 at this epoch), 
"RC" (the epoch at which the surface effective temperature is 10 4 K) and "ffEx" is the epoch at 
which all the free-free photons generated at the epoch marked "ff" are exhausted by radiation from 
the surface. In all cases the epoch is the left most character. Symbols appearing close to either left 
or right vertical axis may have been shifted to keep the symbols within the figure boundary. 



4.2. The Light Curves 

The expected light curves for a macronova at z = 0.2 are shown in Figures 5 and 6. For 
(3 = 0.05 the earliest constraining timescale is "RC" or the epoch marking the recombination 
of the electrons with ions. Beyond day one, the model presented here should not be trusted. 
For P = 0.3, as discussed earlier and graphically summarized in Figure 6, the constraints 
come from photon production. Complete photon exhaustion takes place at the start of epoch 
of transparency (marked by "r es ") and so emission will cease rapidly at 10 4 s. Observations, 
if they have to any value, must be obtained within an hour or two. 

The peak flux in Figures 5 and 6 is about 0.3 /iJy. With 15 minutes of integration time 
on a 10-m telescope one can easily detect this peak flux 

Thus observations are capable of constraining a neutron powered MN with explosion en- 
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Fig. 6. — Expected light curve for a neutron decay powered macronova with j3 s = 0.3. The 
explanations of the symbols can be found in the caption to Figure 5. 

ergy comparable to 10 49 erg (the typical isotropic 7-ray energy release for short hard bursts). 
However, these observations have to be obtained on timescales of hours {(5 S = 0.3) to about 
a day (& = 0.05). 

Observations become less constraining with ejecta speeds larger than that considered 
here because the macronova becomes transparent earlier (t oc assuming that the explo- 
sion energy is constant) and as a result photons are lost on a timescale approaching the light 
crossing timescale i.e. rapid loss. 

A potentially significant and non-trivial complication is confusion of the macronova 
signal by the afterglow emission. Optical afterglow emission has been seen for GRB 050709 
(Hjorth et al. 2005) and GRB 050724 (Berger et al. 2005). Afterglow emission could, in 
principle, be disentangled from the macronova signal by obtaining multi-band data. Of 
particular value are simultaneous observations in the X-ray band. No X-ray emission is 
expected in the macronova model. Thus, X-ray emission is an excellent tracer of the afterglow 
emission. An alternative source for X-rays is some sort of a central long-lived source. As 
discussed in §6 a macronova would reprocess the X-ray emission to lower energies. Thus, 
the X-ray emission, at least in the macronova model, is a unique tracer of genuine afterglow 
emission and as a result can be used to distinguish a genuine macronova signal from ordinary 
afterglow emission. 
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5. Heating by Nickel 

In the previous section I addressed in some detail photon generation and photon-matter 
equilibration for A = Z = 1 plasma. Here, I present light curves with Nickel ejecta. The 
many transitions offered by Nickel (Z = 28, A = 56) should make equilibration less of an 
issue. Next, free- free mechanism becomes increasingly more productive for higher Z ions 
(oc Z 2 ; see Appendix). Thus the time for free-free photons to be exhausted, a critical 
timescale should be longer a Nickel MN (for a given M e j and (3) relative to a neutron MN. 

There is one matter that is important for a Nickel MN and that is the issue of heating of 
matter. In a Nickel MN matter is heated by deposition of gamma-rays released when Nickel 
decays to Cobalt. Relative to the SN case, the energy deposition issue is important for an 
MN both because of a smaller mass of ejecta and also higher expansion speeds. 

The details of Nickel — > Cobalt — > Iron decay chain and the heating function are sum- 
marized in the Appendix. In view of the fact that most of our constraints come from early 
time observations (less than 10 days) we will ignore heating from the decay of Cobalt. Thus 

e m (t) = 3.9 x 10 10 / N iexp(-A N it) ergg'V 1 , (14) 

where A^ 1 = 8.8 d and fai is the mass fraction of radioactive Nickel in the ejecta and is set 
to 1/3. 

Nickel decay results in production of 7-rays with energies between 0.158 and 1.56 MeV 
(see Appendix). Colgate et al. (1980) find that the 7-ray mass absorption opacity is 0.029 cm 2 g 
for either the Ni 56 or Co 56 7-ray decay spectrum. Extraction of energy from the gamma-rays 
requires many scatterings (especially for the sub-MeV gamma-rays) and this "deposition" 
function was first computed by Colgate et al. (1980). However, the e(t) we need in Equation 2 
is the deposition function averaged over the entire mass of the ejecta. To this end, using the 
local deposition function of Colgate et al. (1980), I calculated the bulk deposition fraction, 
?7es (for a uniform density sphere) and expressed as a function of r cs , the center-to-surface 
Thompson optical depth. The net effective heating rate is thus rj es em(t). I find that only 
20% of the 7-ray energy is thermalized for r es = 10. Even for r es = 100 the net effective rate 
is only 70% (because, by construction, the ejecta is assumed to be a homogeneous sphere; 
most of the 7-rays that are emitted in the outer layers escape to the surface). 

The resulting lightcurves are plotted in Figure 7 for f3 s = 0.03 and f3 s = 0.3. I do 
not compute the UV light curve given the high opacity of metals to the UV. As with la 
supernovae, peak optical emission is achieved when radioactive power is fully converted to 
surface radiation (Arnett 1979; Chevalier 1981). 

P s = 0.03 is a sensible lower limit for the ejecta speed; equating the energy released by 



-16- 



radioactive decay to kinetic energy yields (3 S = 0.01. Peak optical emission is achieved at day 
5. At this epoch, the photospheric temperature is in the vicinity of Trc = 5 x 10 3 K. Thus the 
calculation should be reasonably correct up to this epoch. Extrapolating from the neutron 
MN case for similar model parameters I find that photon production and equilibration are 
not significant issues. 

For f3 s = 0.3, peak emission occurs at less than a day, again closely coincident with the 
epoch when the effective photospheric temperature is equal to Trc- For the neutron MN the 
exhaustion of the free-free photons was the limiting timescale (~ 10 4 s). Qualitatively this 
timescale is expected to scale as Z (§5) and if so the model light curve up to the peak value 
is reliable. 




Fig. 7. — Model light curve for a macronova at z = 0.22 and powered by radioactive Nickel (one 
third by mass). The dotted line is the expected light curve if Nickel radioactive decay power is 
instantly converted to the optical band. The input includes only the fraction of 7-rays that are 
absorbed and thermalized within the ejecta (i.e. £Ni(^) r /es)- The slight curvature in the dotted line 
(between 10 2 s and 10 4 s) is an artifact of the least squares fit to the net energy deposition function 
(Equation 47). "RC" is the epoch at which the effective surface temperature is 5,000 K and the 
epoch at which the electron optical depth is 10 is marked. 

To conclude, a Nickel powered MN is detectable only if the explosion speed is unrea- 
sonably low, f3 s = 0.03. Observations will not place significant constraints for f3 s = 0.3. 
For such rapid expansion, the MN suffers from lower deposition efficiency and an onset of 
transparency before the bulk of Cobalt has decayed. Conversely, there exists an opportunity 
for (futuristic) hard X-ray missions to directly detect the gamma-ray lines following a short 
hard burst! 
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6. The MN as a reprocessor 

The simplest idea for short hard bursts is a catastrophic explosion with an engine that 
lives for a fraction of a second. However, we need to keep an open mind about the possibility 
that the event may not necessarily be catastrophic (e.g. formation of a millisecond magnetar) 
or that the engine or the accretion disk may live for a time scale well beyond a fraction of a 
second. 

It appears that the X-ray data for GRB 050709 are already suggestive of a long lived 
source. A strong flare flare lasting 0.1 d and radiating 10 45 erg energy (and argued not to 
arise in the afterglow and hence by elimination to arise from a central source) is seen sixteen 
days after the event (Fox et al. 2005). The existence of such a flare in the X-ray band limits 
the Thompson scattering depth to be less than unity at this epoch. As can be seen from 
Figures 5 and 6 this argument provides a (not-so interesting) lower limit to the ejecta speed. 

The main value of a macronova comes from the reprocessing of any emission from a 
long-lived central source into longer wavelengths. In effect, late time optical observations 
can potentially constrain the heating term in Equation 2, regardless of whether the heating 
arises from radioactivity or a long-lived X-ray source. In this case, e(t)M e j refers to the 
luminosity of the central source. The optical band is the favored band for the detection of 
such reprocessed emission (given the current sensitivity of facilities). 

A central magnetar is a specific example of a long lived central source. The spin down 
power of an orthogonal rotator is 

where B is the dipole field strength, R n is the radius of the neutron, uj = 2n/P is the rotation 
angular frequency and P is the rotation period. For B = 10 15 G, R n = 16 km we obtain 
dE/dt ~ 10 42 (P/100 ms)~ 4 erg s" 1 and the characteristic age is 5 x 10 4 s. Constraining such 
a beast (or something similar) is within the reach of current facilities (Figure 8). 

The power law decay discussed in this section is also applicable for neutron rich ejecta. 
For such ejecta arguments have been advanced that the resulting radioactive heating (from 
a variety of isotopes) can be approximated as a power law (see Li & Paczyhski (1998)). 

Earlier we discussed of the potential confusion between a macronova and afterglow 
emission. One clear distinction is that afterglow emission does not suffer from an exponential 
cutoff whereas a macronova dies down exponentially (when the optical depth decreases). 
However, quality measurements are needed to distinguish power law light curves of afterglow 
(flux oc t a with a — — 1 to —2) versus the late time exponential cutoff of macronova emission. 
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Fig. 8. — Expected optical light curve from a power law central source with L m (t) = L m (0)(l + 
t/tmiO) 2 )' 1 with L m (0) = 2 x 10 41 erg s _1 and t m (0) = 10 5 s. The dotted line is the expected flux 
if the power from the central source is instantly and fully converted to the optical band. "RC" and 
"r = 10" have the same meanings as in previous figures. 



The prevailing opinion is that short duration bursts arise from the coalescence of a 
neutron star with another neutron star or a black hole. The burst of 7-rays requires highly 
relativistic ejecta. The issue that is very open is whether the bursts are also accompanied 
by a sub-relativistic explosion. This expectation is motivated from numerical simulations 
of coalescence as well as the finding of supernovae accompanying long duration bursts. Li 
& Paczyhski (1998) were the first ones to consider the detection of optical signal from sub- 
relativistic explosions. Rather than referring to any such accompanying nova-like explosion as 
a "mini-supernova" (a term which only an astronomer will fail to recognize as an oxymoron) 
I use the term macronova (and abbreviated as MN). 

Essentially a macronova is similar to a supernova but with smaller mass of ejecta. 
However, there are important differences. First, supernova expand at relatively slow speed, 
~ 10 9 cms _1 (thanks to the smalle escape velocity of a white dwarf). Next, radioactivity, 
specifically the decay of radioactive Nickel plays a key role in powering supernovae. The 
large mass of the ejecta, the slow expansion and the 9-day and 111-day e-folding timescale of 
Nickel and Cobalt (resulting in a gradual release of radioactive energy) all conspire to make 
supernovae attain significant brightness on a timescale well suited to observations, namely 
weeks. 



7. 



Conclusions 
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A macronova has none of these advantages. On general grounds, the mass of ejecta is 
expected to be small, 1O~ 2 M . The expansion speed is expected to be comparable to the 
escape velocity of a neutron star, 0.3c. Finally, there is no reason to believe that the ejecta 
contains radioactive elements that decay on timescales of days. 

In this paper, I model the optical light curve for a macronova powered by decaying 
neutrons or by decay of radioactive Nickel. The smaller mass and the expected higher 
expansion velocities necessitate careful attention to various timescales (photon generation, 
photon-matter equilibrium, gamma-ray energy deposition). 

Surprisingly a neutron powered MN (with reasonable explosion parameters) is within 
reach of current facilities. Disappointingly a Nickel powered MN (with reasonably fast ex- 
pansion speed) require deeper observations to provide reasonable constraints. This result 
is understandable in that a MN is not only a smaller supernova but also a speeded up su- 
pernova. Thus the detectability of a MN is intimately connected with the decay time of 
radioactive elements in the ejecta. Too short a decay (seconds) or too long a decay (weeks) 
will not result in a bright signal. For example, artificially changing the neutron decay time 
to 90 s (whilst keeping the energy released to be the same as neutron decay) significantly re- 
duces the signal from a macronova so as to (effectively) render it undetectable. The difficulty 
of detecting Nickel powered MN illustrate the problem with long decay timescales. 

Next, I point out that a central source which lives beyond the duration of the gamma-ray 
burst acts in much the same way as radioactive heating. 

Finally, it is widely advertised that gravitational wave interferometers will provide the 
ultimate view of the collapse which drive short hard bursts. However, these interferometers 
(with appropriate sensitivity) are in the distant future. Rapid response with large telescopes 
has the potential to directly observe the debris of these collapses - and these observations can 
be done with currently available facilities. This scientific prize is the strongest motivation 
to mount ambitious campagins to detect and study macronovae. 

To date, GRB 050509b was observed rapidly (timescales of hours to days) with the most 
sensitive facilities (Keck, Subaru and HST). The model developed here has been applied to 
these data and the results reported elsewhere. 

This paper is the first theoretical (even if mildly so) paper written by the author. The 
pedagogical tone of the paper mainly reflect the attempts by the author (an observer) to make 
sure that he understood the theoretical underpinnings. I gratefully acknowledge pedagogical 
and fruitful discussions with P. Goldreich, R. Sari and S. Sazonov. I would like to thank 
L. Bildsten for acting as the internal referee, U. Nakar for help in understanding the ODE 
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solver in MATLAB, R. Chevalier, T. Piran, A. MacFadyen, B. Schmidt and R. Sunyaev for 
suggestions and clarifying discussions. The author is grateful for Biermann Lecture program 
of the Max Planck Institute for Astrophysics, Garching for supporting a one month sabbatical 
stay. I am very grateful to T.-H Janka for patient hearing and encouraging me to submit this 
paper. The author acknowledge financial support from a Space Telescope Science Institute 
grant (HST-GO-10119) and the National Science Foundation. 
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Appendix: Formulae & Approximations used in the paper 

Blackbody Radiation 

The blackbody spectral intensity is 

B "< T > = i eM-hWkT) - 1 ^cm-V'ster-'Hz-'. (16) 
Frequently it helps to use physically normalized frequency unit, x = hvj{k B T). Then, 



The spectral energy density is 



The photon number density is 



u u {T) = —B V (T). 
c 



„«„) = / (18) 

c hv 



n[ < x) 



= 8n(^Y [ X * dx (19) 
\ h J J exp(x) — 1 

= Sn(^) 3 m< 3 (x) (20) 

= ) 3 C3(,); (21) 



This preprint was prepared with the AAS IATgX macros v5.0. 
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where ( n (x) = r(n) -1 J* t n_1 /(e* — l)dt is an incomplete ( function of order n. We note 
that C3(°°) = 1-202 and C4(°°) = tt 4 /90. Thus the number of photons per unit volume is 
n = 19.2n(kT/hc) 3 . The mean energy of the photon defined as u(T)/n is 2.7 kT. In the 
Rayleigh- Jeans approximation, we obtain 

„« ,) = 16 ,(_) _ (22) 

This equation demonstrates that most of the photons have x > 1 (with less than 20% of 
photons with x < 1). The fiftieth percentile is x = 2.36. 



Free-free Emission 

The free-free volume emissivity rate (Rybicki & Lightman 1979, p. 160) is 

€b(u) = 6.8 x 10~ 38 Z 2 n e niT~ 1/2 exp(-hv/k B T)g s erg cm^s^Hz" 1 . (23) 

Here, g& is the velocity-averaged Gaunt factor and varies from 5 to 1 as u = hu/k B T ranges 
from 10~ 4 to 1. It appears to me that = — log(w) + 1 provides a reasonable approximation 
(for u < 1). The frequency integrated emissivity (Z — A — 1) 

e ff = 1.4 x 10~ 27 n 2 e T l/2 g B erg cm" 3 s" 1 (24) 

where g B is the frequency-averaged Gaunt factor and a value of 1.2 gives 20% accuracy. For 
very hot plasma, relativistic effects become important and the rate increases by a factor of 
1 + T/2.2 x 10 9 K. 

The free-free density absorption coefficient, a u (&) = na l/ (S) is given (it op cit. , p. 162) 
a u (S) = 3.7 x 10 8 T- 1/2 Z 2 n e n^~ 3 (l - exp(-hv/k B T)) cm" 1 ^ (25) 

As noted in the previous section, it is useful to work in normalized frequency units, 
x = hvjkT. Noting e s (x) = {kT/h)e s {u) we find (Z = A = 1) 

e s (x) = 1.42 x 10~ 27 n 2 T 1/2 exp(-ar) erg cm' 3 s" 1 . (26) 

The number of photons above a value of x is given by the integral y(x) = t" 1 exp(— t)dt. 
Though simple looking, this integral has no solution with known functions. The usual (sen- 
sible) approximation recognizes that there are few photons with x > 1 and thus the upper 
limit can be replaced by unity. The next approximation is replace the exponential by unity. 
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Thus y(x) = —\og(x) (with implicit understanding x < 1). Alternatively, an empirical 
function (obtained by numerical fitting), 

y(x) = 0.306 - 0.413 log(x) + 0.094 log(:r) 2 , (27) 

provides a reasonable fit for 10~ 2 < x < 5. 



Photon-electron equilibration timescale 

For a pure A = Z = 1 gas, the primary primary photon-matter interaction for very 
hot gas is through Thompson scattering of electrons. Here we summarize the basic physics 
of the timescale of thermalizing a fast moving electron by the radiation field as well as the 
thermal photon-thermal electron equilibration timescale. 

Consider the fate of the electron ejected as a result of neutron decay. The initial Lorentz 
factor of the electron (70) can range from 3.5 to unity with a mean value of 1.6. The typical 
ejected electron is much hotter than the radiation field (temperature T) and will cool down 
by inverse scattering the ambient photons. The rate of energy loss is 

d -§ = (28) 

where the kinetic energy of the electron at any instant is given by E = (7 — l)m e c 2 . Simpli- 
fying I obtain 

<N = A(l - 7 2 ) (29) 



dt 

where A = AarcU-^/ZmeC 2 . Since J dx(l — x 2 )" 1 = (1/2) log((x + l)/(x — 1)) we find that 
the electron loses a large fraction of energy over timescales of A -1 . 

Next we consider an electron with energies comparable to that of photons. Consider 
an electron moving with velocity v with respect to the radiation field. In the electron rest 
frame the radiation temperature is Doppler shifted as follows: 

T(0) =T(l + v/c)cos(6) (30) 

and correspondingly the total intensity also has an angular dependence given by 

1(0) = ^-T{6)\ (31) 

The momentum is 1(9) jc and the rate of momentum transfer to the electron is 

F = - J <T r p—^—- cos(9)dQ (32) 
= -\<y T aT'- (33) 



-25 - 



Averaging the velocity over the thermal distribution, the work done by the radiation 
drag on the electron is 

(-Fv) = ^ T T 4 (^) (34) 
kT 

= 4aTaT 4^£ (35) 
m e c 

where we have made use of the well known formula l/2m e (v 2 ) = 3/2kT e . However, the elec- 
trons must be getting heated up by fluctuations in the radiation field because in equilibrium, 
when T = T e , the electron will neither gain energy nor lose (on average). Thus the heat 
equation for the electron is 

^ = 4a T aT 4 k(T - T e )/m e c. (36) 

For fully ionized matter, Q = 3kT e and thus we get the equation well known to cosmologists 
(Peebles 1971) 

dT 

= -A(T e - T) (37) 
where U p h = aT 4 . The thermalization timescale is thus ~ A^ 1 . 



Electron and Ion Equilibration 

A fast electron will lose energy to slow moving electrons by Rutherford scattering. From 
Padmanabhan (2000, volume I, pp. 439) we see 

«*•*-> - < 38 > 

t(i,i) = — : = 4.3x10 T i 7 n u s, and (39) 

^fc/T fc- f vj-ji^ 

i(e ' !) = {T)' /2 H^y M=wz ~ 2 ^ n " s - (4o » 

Here M is the mass of the ion and L e (Lj) refers to the usual logarithmic factor for electrons 
(ions) and set to 10. n is the particle density (cm -3 ) and the density normalization of 
10 11 cm -3 is appropriate for radius of 10 14 cm and ejecta mass of M cj = 10~ 3 M Q . 



Heating by Neutron Decay 



A free neutron decays to a proton, electron and antineutrino. The half lifetime of the 
decay is 10.4 minutes and thus the e-folding time, t N = 900 s. The mass difference between 



-26- 



the neutron and that of the electron and proton is Q = 0.782 MeV (see for example, Krane 
1988). 

The sum total linear momentum is conserved and the energy carried off by the electron 
and antineutrino must add up to Q. This would nominally suggest that there are five 
remaining free parameters (3x3 momenta minus four constraints). Thus the output of decay 
can now be described by five free parameters and we choose these to be the momentum of 
the electron, p, and two momenta of the antineutrino, q. 

The proton, by the virtue of its mass, obtains very little energy in the process. Thue the 
energy of the antineutrino is E p = Q — T where, following the convention in nuclear physics, 
T = E — m e c 2 is the kinetic energy and E 2 = p 2 c 2 + m 2 c A (here p = y/p ■ p) . The magnitude 
of the momentum of the antineutrino is q = E p /c (since the antineutrino has very low, if 
any, mass) 

There is no preferred direction for p and the phase space (per unit volume) for the 
electron momentum to lie between p and p + dp and the antineutrino momentum to lie 
between q and q + dq is 

4irp 2 Anq 2 

dn e dn 9 = —r^-dp—-—dq. (41) 
h 6 h 6 

Our goal is to determine the distribution of the momentum of the electron. To this end we 
integrate dn e drip over q but subject to the constraint qc = Q — T and find 

dn e (p) otp 2 (Q -T) 2 dp. (42) 

We transform this distribution to the differential energy distribution of the electron, dn e (E)/dE, 
using the usual rule of transformations 

dn e (E) dn e (p) dp 

dE ~ dp ~dE 

oc Ep(Q-T) 2 . (43) 

For the case of a decaying neutron, we use the probability distribution given by Equa- 
tion 43 and find that the mean kinetic energy of the electron is (T) = 0.59 m e c 2 . Thus the 
heating rate per gram is Na{T) /t^ = 3.2 x 10 14 ergg _1 s -1 . 



Heating by Decay of Radioactive Nickel 

56 Ni is produced whenever nuclear reactions are allowed to proceed to statistical equi- 
librium. This isotope enjoys the distinction of being the most packed nucleon (highest 
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binding energy per nucleoli). However, the isotope is unstable and decays as follows: 
56 Ni(6.01 d)— » 56 Co(77d)— » 56 Fe. The half lifetimes are given in the parenthesis. The de- 
cay proceeds by electron capture of a K shell electron. Following this capture, a neutrino is 
emitted and the 56 Co is left in an excited state, 1.72 MeV above the ground state. A variety 
of decays are possible from this state to the ground state of 56 Co: 1.56 MeV, 0.812 MeV, 
0.75 MeV, 0.48 MeV, 0.27MeV and 0.158 MeV (see Arnett 1996, p. 426). Cobalt decays by 
electron capture (80%) resulting in a spectrum of 7-rays of mean energy of 3.5 MeV per 
decay and beta-emission (20%) with mean energy of 0.14 MeV per decay (see Colgate et al. 
(1980)). We set (E Co ) = 3.5 + 0.14 = 3.64 MeV. 

Let N = N-Ni{t = 0) be the initial number of Nickel nuclei. Let Nc be the number of 
Cobalt nuclei; we have iV Co (t = 0) = 0. The time evolution of Nickel and Cobalt is given by 
well known parent-daughter relations: 

Nm(t) = N exp(— \mt) and 

N Co (t) = N Xm fexp(-A C ot)-exp(-A Ni t)Y (44) 

^Ni — ^Co V / 

Here, Ani = t^l and Ac = ^co with t^i = 6.1/ log(2) d and tc = 77/ log(2) d. The heating 
rates (per gram) are obtained by letting N = (56m#) _1 and are 

£Ni(f) = XmAE m N Ni (t) 

= 3.9 x 10 10 exp(-A Ni t) ergg^s" 1 

s C o(t) = Xco(Eco}N C o(t) 

= 7 x 10 9 ^exp(-A Co t) -exp(-A Ni t)^ ergg^s" 1 . (45) 

From an inspection of the heating curve one finds that Cobalt decay competes with Nickel 
(assuming full energy absorption) by day 17. In the main text we note that early time (<; 1 
week) observations result in the best constraint (or even a discovery). Given this we can 
safely ignore the heating provided by the decay of Cobalt. 

The next issue is efficiency with which the 7-rays emitted during the course of the decay 
of Nickel can transfer energy to the surrounding matter. Colgate et al. (1980) consider the 
deposition of energy by 7-rays emanating from decay of Ni 56 and Co 56 . They find that the 
mass opacity of the 7-ray absorption is about 0.029 cm 2 g -1 , for either Ni 56 or Co 56 decay 
spectrum. They used a Monte Carlo photon transport code and found that at least for 
a uniform sphere, D, the fraction of the 7-ray absorbed by the matter at a given radius 
depends only on 7-ray optical depth (to the surface), r 7 : 



£>(r 7 ) = G[l + 2G(1 - G)(l - 0.75G)] 



(46) 
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where G = r 7 /(1.6 + r 7 ). 

I integrated the net deposition fraction for spheres with center-to-surface optical Thomp- 
son optical depth, r es ranging from 1 to 10 17 . A polynomial fit yields 

logics)) = -0.00032 4 + O.OIO82 3 - 0.1486z 2 + 0.8134,2 - 1.4360 (47) 

where z = log(r es ); here log refers to logarithm to the base 10. 

For example, when r es = 10 I find 7] 1 = 0.2 i.e. 80% of the 7-rays escape to the surface. 
Even for r es = 100 I find r? 7 = 0.7 i.e. about 30% escape to the surface (mainly the outer 
layers) . 



